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We analyze the interaction of plane '+'-polarized gravitational waves with a Fabry-Perot cavity in 
the local Lorentz frame of the cavity input mirror outside of the range of long-wave approximation 
with the force of radiation pressure taken into account. The obtained detector response signal is 
represented as a sum of two parts: (i) the phase shift due to displacement of the movable mirror 
under the influence of gravitational wave and the force of light pressure, and (ii) the phase shift 
due to direct interaction of gravitational wave with light wave inside the cavity. We obtain formula 
for the movable mirror law of motion paying close attention to the phenomena of optical rigidity, 
radiative friction and direct coupling of gravitational wave to light wave. Some issues concerning the 
detection of high-frequency gravitational waves and the role of optical rigidity in it are discussed. 
We also examine in detail special cases of optical resonance and small detuning from it and compare 
our results with the known ones. 
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I. INTRODUCTION 

The search for gravitational waves (GWs) is now con- 
ducted with large Earth-based Michelson-type laser inter- 
ferometers where test masses are separated by a distance 
of several kilometers [H, 0, 9 ■ Variations in the proper 
distances between these test masses are measured with 
the light wave. Traditionally interferometric GW detec- 
tors are analyzed in the transverse-traceless (TT) gauge, 
where coordinate lines are chosen to coincide with the 
geodesies of the freely falling test masses (if initial veloc- 
ities of test masses are equal to zero) [3, [fl ■ This means 
that the coordinate frame is rigidly bounded with the in- 
ertial test masses. For a '+'-polarized plane gravitational 
wave metric in the TT gauge has the usual linearized 
form 

ds 2 —g a /}(x)dx a dx /3 = [r) a p + h a p(x)] dx a dx@ 

= - c 2 dt 2 + (1 + h)dx 2 + (1 - h)dy 2 + dz 1 . (1) 

Here h — hit — z/c), \h\ <C 1 is the GW wavefunction, 
Va0 = diag(— 1, 1, 1, 1) is the metric tensor of flat space- 
time, and h a p are small deviations from it. Greek indices 
run over 0, 1, 2, 3 or ct, x, y, z\ Latin indices run over 1, 2, 3 
or x, y, z. 

Geodesic equation in metric dXJ) dictates that the test 
mass coordinates are not changed by gravitational wave 
[j|, i.e. GW does not interact with the test masses 
in the TT gauge. Nevertheless, it noticeably changes 
the phase of electromagnetic wave (EMW for short) 
which propagates in the region of space-time where grav- 
itational wave is present. The elementary solution of 
eikonal equation for a single arm of laser interferome- 
ter Q states that the phase shift is proportional to GW 



signal: <5\E'(t) = k^Lhit), where kg is the wavenumber 
of electromagnetic wave and L is the interferometer arm 
length. In other words one can reduce the analysis of a 
'+ '-polarized GW interaction with detector to the prob- 
lem of EMW propagation in the medium with effective 
refraction index n 2 (t) = 1 + h(t). This result is only 
valid in long- wave approximation L <C A gw , where A gw 
is the length of gravitational wave. In general case the 
formula for phase shift is more complex but it is still the 
solution of the eikonal equation. The significant advan- 
tage of the TT gauge is the ease of going from long-wave 
approximation to the general case [5(. 

However, the crucial requirement of the test masses 
inertiality cannot be met in the real Earth-based ex- 
periments since there exists a great amount of non- 
gravitational forces acting on the test masses, such as 
seismic and thermal fluctuating forces, light pressure re- 
sulting in pon deromotive rigidity @, [WE, III [II [3 
EH [3 EH H3 an d radiative friction [18l]T The latter two 
phenomena appear if Fabry-Perot (FP) cavities are uti- 
lized in the interferometer as in LIGO and VIRGO detec- 
tors 3 . Therefore, the reference frame of the TT gauge 
can hardly be constructed for the Earth-based detectors. 
Correspondingly, the results obtained within the bounds 
of the TT gauge are only applicable when condition of 
the test masses inertiality is met (the case of space GW 
antenna LISA [lj|, [2(| in which drag-free technique will 
be used). Otherwise, additional analysis should be per- 
formed in order to validate the TT-based consideration 
of non-geodesic motion. Clearly, this will significantly 
complicate the whole problem. 

In order to get rid of inertiality requirement when an- 
alyzing the interaction of GW with a Fabry-Perot cavity 
and take non-gravitational forces into account we should 
work in the local Lorentz frame (LL gauge) of the cavity 
input mirror. Coordinate transformation from the TT 
gauge to the LL gauge can be found in Ref. [f|: 
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ds 2 = g a p(x)dx a dx (3 = - c 2 dt 2 + dx 2 + dy 2 + dz 2 
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+ ^ (x 2 -y 2 )h(cdt- dzf, 

(2) 

where old notations for the coordinates are used since 
we are not interested in the coordinates of the TT gauge 
anymore. 

A simple analysis performed in Ref. [j| shows that in 
the LL gauge gravitational wave mainly interacts with 
the test masses, resulting in the phase shift S^(t) = 
2ko5X(t) = koLhif), while direct interaction with the 
light wave leads only to a small correction in long-wave 
approximation. Here SX(t) = (l/2)Lh(t) is the displace- 
ment of the test mass due to its interaction with gravita- 
tional wave. In other words one can reduce the analysis 
of a '+ '-polarized GW interaction with detector to the 
problem of EMW propagation in the volume without ef- 
fective medium (in long-wave approximation only!) but 
with movable boundaries. 

However, the LL gauge has one great disadvantage — 
the complexity of going from long- wave approximation to 
the general case. In this paper we are mostly interested 
in the latter because still little is known about the partic- 
ular features of non-inertial test masses motion under the 
influence of gravitational wave when condition L <C A gw 
is violated. In this sense we are mostly interested in the 
phenomenon of optical rigidity (or optical spring) and 
its role in detection of high frequency GW signals. It is 
most likely that the optical strings will be implemented in 
the future generations of the Earth-based detectors (such 
as Advanced LIGO), allowing significant amplification of 
the response to GW signals. General formulas for optical 
rigidity, test masses motion and detector response signal 
outside of the range of long- wave approximation were ob- 
tained in Ref. [ll| . Starting from the ideas of Ref. [2l| 
we propose an independent derivation of general formu- 
las and obtain a good agreement with results of Ref. [ll| 
on the one hand, and of [HI, [l8j] on the other. Our cal- 
culations allow us to take into account the relativistic 
correction to the force of radiation pressure, known as ra- 
diative friction [l8j] , which was not obtained in Ref. [ll| . 
Using modern notations we further analyze the GW- and 
radiation pressure-driven motion of test masses, paying 
close attention to general relativistic aspects of the prob- 
lem. We also study in detail the differences between the 
general case and its long-wave approximation. It should 
be noted that we deal only with the optimal GW polar- 
ization and detector orientation ('-(-'-polarized GW with 
one of its principal axes parallel to the axis of FP cav- 
ity) since it is enough for consideration of the radiation 
pressure-induced dynamics. 

Motivation for our work is the following. Modern 
GW detectors, utilizing Fabry-Perot cavities, with arm 
lengths of several kilometers (4 km for LIGO and 3 km for 
VIRGO) are best suited for detection of GW signals with 
frequencies that lie in the region of 10 2 to 10 3 Hz. For 
them L/Ag W ~1CT 3 -^10~ 2 <S;1 and long- wave approxi- 
mation holds with high accuracy. Today the most efforts 
are made to detect these low-frequency signals which are 



emitted, for example, by binary pulsars, rotating neutron 
stars, supernova explosions etc. However, modern astro- 
physical and cosmological theories predict the existence 
of sources of high-frequency (/ > 10 4 Hz) GW signals, 
such as relic gravitational wave background [13, HH, [24[ • 
Some models, such as inflation, predict the GW spectrum 
of the background with typical values of amplitude which 
are far beyond the sensitivity of even advanced interfer- 
ometers. However, the spectrum of g ravitational waves 
predicted by string cosmology [HI, [26[[27| has rather big 
typical GW amplitude which rises with increasing fre- 
quency, thus allowing the potential possibility of detec- 
tion of GWs with certain frequencies above 10 3 Hz on 
the Earth-based detectors. There is a strong evidence 
that modern and future detectors are (or will be) capa- 
ble of detecting signals at frequencies about 10 4 Hz. For 
example, in Refs. [28L [29l| (see also Refs. [3(1 HU) ^ was 
proposed that 4 km LIGO detectors could detect 37.5 
kHz signals, corresponding to cavity free spectral range 
(FSR), for which L/A gw = 1/2, thus operating outside of 
long-wave regime. And though in this case the response 
signal vanishes for optimal polarization of gravitational 
wave [Hj] , GW signals can still be detected with enhanced 
sensitivity at multiples of the cavity FSR for other ori- 
entations. However, the role of radiation pressure effects 
in the detection of high-frequency GWs when condition 
L I A gw <C 1 is violated (especially in the vicinity of FSR 
and its multiples) has not been recognized up to date. 
We propose the study of these effects, since the searches 
for GW background are planned for the Advanced LIGO 
stage where optical springs will be utilized. 

This paper is organized as follows. In Sec. [II] we write 
down the wave equation for electromagnetic wave propa- 
gating in the space-time of plane gravitational wave and 
obtain its solution for an arbitrary h(t). In Sec. IIIII we 
examine the special case of monochromatic h{t). In Sec. 
HVl we apply these results to consideration of the bound- 
ary problem for EMW in a Fabry-Perot cavity with the 
movable mirror and obtain its response signal to GWs. 
By neglecting the resonant features of FP cavity in Sec. 
[V]we find excellent agreement between our LL-based cal- 
culations and TT-based calculations performed in Ref. 
Q for the case of single reflection. In Sec. IVII we write 
down full equation of motion for the mirror paying close 
attention to the phenomena of optical rigidity, radiative 
friction and GR corrections to the force of radiation pres- 
sure. In Sec. IVIII we find the mirror law of motion and 
substitute it into the response signal. The obtained for- 
mulas basically agree with results of Ref. [ll| . We also 
discuss the possibility of parametric excitation of addi- 
tional optical modes under the influence of GW signals 
with frequencies near cavity FSR and its multiples. In 
Sec. IVIIII we examine special cases of optical resonance 
and small detuning from it in general case and long-wave 
approximation. Excellent agreement with the results ob- 
tained in Refs. [H, [Hj] is found. In the end we estimate 
the deposit of optical rigidity into the response signal in 
the vicinity of FSR. 
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II. WAVE EQUATION AND ITS SOLUTION 

We start from the second pair of Maxwell's equations 
without sources in curved space-time 



1 



d 



-g dx 13 



a/3\ 



0. 



Here F MV = d fJ ,A v - d v A^ A» = (A , A\ A 2 , A 3 ) is the 
4-potential of electromagnetic field, and g = det(g PlL ,) = 
— 1. Let us impose the Coulomb gauge A = 0, diA 1 + 
d 2 A 2 + d 3 A 3 = 0. Then 

g a0 d a d A^ = 0. 

In general A^ = A p, (x,y, z,t). Let us restrict ourselves 
to consideration of gravitational waves with frequen- 
cies / gw ~ 10 2 -j- 10 5 Hz, corresponding to wavelengths 
A gw ~ 10 3 -j- 10 6 m. This is far larger than the opti- 
cal wavelength Ao ~ 1 mkm and the radius of the laser 
beam r ss 10 cm. For simplicity we assume that A^ de- 
scribes a linearly z-polarized plane light wave traveling 
along the £-axis. Therefore, under all listed conditions 
we have A^ = (0,0,0, ,4) with A = A(x,t) and the fol- 
lowing scalar wave equation: 



d 2 A 1 d 2 A 



dx 2 



c 2 dt 2 



1 x 2 •■ , , 1 d 2 A 
h{t) 



2 c 2 



c 2 dt 2 



(3) 



The right side of this equation describes the direct inter- 
action of gravitational wave with light wave [5(. In Ref. 
[2~fl ] this effect was called the distributed gravitational 
redshift. 

It is convenient to solve Eq. using the method 
of successive approximations since h <C 1. We shall 
keep only the zeroth and the first order in h terms: 
A(x,t) = A(°\x,t)+AW(x,t), ~ \hAW\ < 

The solution of the zeroth order is the sum of plane 
monochromatic waves traveling in positive and negative 
directions of the x-axis. We denote "positive" wave with 
index '+' and "negative" wave with index '-': 

A^(x,t) = A { °\x,t)+A {0) (x,t), 

A ( i ] (x, t) = A ±oe - l{uJotTkox) + c.c, (4) 

where fco — wo/c. Amplitude and frequency are derived 
from some initial and boundary problems and we shall 
keep them undefined until Sec. HVl The first order equa- 
tion is: 



8 2 A^ _ 1 d 2 A^ _lx 2 .. 1 d 2 A^ 
dx 2 ~ c 2 dt 2 ~ 2 c 



2^)-^^- (5) 



z 2 dt 2 ' 

The general solution of this equation can be represented 
as the sum of "positive" and "negative" waves: 



A« (x, t) = A { 1> {x, t) + A w (x, t) 



(6) 



Clearly, they can be treated independently. Remind, that 
goo(0,t) = — 1 and therefore, we must demand that 



AV>(p,t)=AW(Q,t) = 0. 



(7) 



Physically these initial conditions mean that both the 
light waves A±(x, t) experience no redshift at x = 0, 
i.e. the solution of full Eq. ^ at x = is A±(0,t) = 

A^\0,t) = A± Q e-' lulat + A* ±0 e luJot . The solution of the 
Cauchy problem (J5HZJ) is obtained in Appendix 



A±\x, t) = A ±0 g ± (x,t)e^ ot ^ + 



c.c, 



(8) 



+ 00 



g±(x,Cl + w )e 



-int 



2^' 



where 



g±(x, Q) = - h(Q - ujq) 

x \ x 2 £_(fl) =F ixr](Q) - ((ft) 



and 



L0 



7,(0)^4^ 



c 2 f2 + luq 

LU 



C(o) 



m 2 (fl 2 + 3lu 2 ) 



(9a) 
(9b) 
(9c) 



III. 



CASE OF MONOCHROMATIC 
GRAVITATIONAL WAVE 



In order to make clear the physical meaning of the 
obtained result we shall examine the special case of 
monochromatic gravitational wave h(£i) — 2nho5(fl — 
Wgw) + 2irhQ5(fl + w gw ), corresponding to h(t) = 
2ho cos u>g W t in time domain. We shall restrict our- 
selves to consideration of the light wave traveling in pos- 
itive direction of the x-axis and use natural requirement 
w gw /27r ~ 10 2 H-10 5 Hz < w /2tt - 10 14 Hz. Under these 
conditions formulas (|4]) and |8|) lead to: 

A+{x,t) = A+{x,t)e l{koX - Uot) ^ +{x ^ 



= A, 



1 — ihr 



a(x) sinWg W t + b(x) cosuj gw t 



('-!) 



where 



1 x 



CJ 



a ( x ) = 77 ^ w 0^gw 

2 & to™ 



>e i(k x-uj t) , c c 



b{x) = ~Uq. 
c 



We can now easily extract information about the ampli- 
tude and the phase of light wave. Up to the second order 
of h amplitude remains unchanged A+(x, t) — A + q while 
phase "i+ix, t) is influenced by GW in the first order: 



*+(ar,t) 



• LJqLO 



gw 



sin tu gw t 
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X Wo f x\ 
— u)q cos Lu gV! t H sm w gw I r J 



IV. LIGHT WAVE IN A FABRY-PEROT 
CAVITY UNDER THE INFLUENCE OF PLANE 
GRAVITATIONAL WAVE 



We must put t = x/c if we are interested in the phase of 
light wave at point x. This is enough since taking into 
account the influence of gravitational wave on time will 
lead to the term of the second order in h. 



^> + (x,x/c) = -h Q 



X UJq 
CJO^gw 



X X 
— UJq COS UW — 

c c 



We plotted the graph of this function in Fig. |T|). 



x 10" 




FIG. 1: Phase shift ^ +(x,x / c) acquired by "positive" wave 
from its direct interaction with monochromatic gravitational 
wave. For visualization we assumed u>o/2n = 3 x 10 14 Hz, 
u gv /2ir = 37.5 kHz, ho = 10" 22 . 



The long- wave approximation x / A gw C 1 of this for- 
mula is 



*+(*) 



— - — k xh -rj- + k xh — — 



O knxho 



A 2 

*gw 



(10) 



This result coincides with the one of Ref. [5[: the di- 
rect interaction of gravitational wave and light wave leads 
only to a small correction of the order of (x/A gw ) 2 h in the 
LL gauge. Note that formula (fTU)) can be qualitatively 
interpreted as the propagation of light wave in effective 
medium with refraction index n 2 (x,t) = 1 + 5n 2 (x,t), 
where Sn(x, t) cx (x/c) 2 h(t), x <C A gw . 



Let us consider now the light wave traveling inside the 
Fabry-Perot cavity placed in the space-time of a plane 
'+ '-polarized gravitational wave. The cavity is formed 
by two mirrors (see Fig. The first one (left) is called 
the input mirror; it is immovable in its own reference 
frame (remind that we work in the local Lorentz frame 
of this mirror) and partially transmits optical radiation 
with amplitude coefficient T< 1. The second mirror — 
the end-mirror — is movable and absolutely reflective. 
For simplicity we assume that both mirrors are ideal, 
meaning that they do not absorb or dissipate optical en- 
ergy. We put distance between the mirrors in the absence 
of gravitational wave and optical radiation to be equal to 
L. 



A[ n (x,t) 



R, T 



^OUtf^) t) 



A + (x,t) 



LA 



A.(x,t) 



R=l 

A \ i 
1 1 
1 1 

ZJ L.i 



L L+X(t) 

FIG. 2: Fabry-Perot cavity with the movable mirror. Input 
wave is denoted by index 'in', output wave is 'out'. Optical 
field inside the cavity is a sum of two waves running in two 
opposite directions of the :r-axis: "positive" wave is denoted 
by index '+' and "negative" is '-'. Isolated cavity without any 
gravitational or optical radiation present has a fixed length L. 
Coordinate of the movable mirror is L + X(t), Ici. 



Let the plane linearly polarized light wave j4i n (a;, t) be 
incident on the input mirror. We shall write it down as a 
sum of "strong" wave with amplitude Aj n o and frequency 
loq and "weak" addition a- ln (x,t) which corresponds to 
the noises of light field: 

A- m {x, t) = A- mQ e-^ ot - koX *> + a in (x, i) e -K-ot-/c„*) + c ^ 
a in {x,t) = f a in (r! + wo)e- ia ( t ~ f ) ^. 

For simplicity in this paper we consider classical EM 
fields (and corresponding classical noises), since we are 
focused on dynamics only. The generalization for quan- 
tized EM fields is straightforward: classical field ampli- 
tudes should be replaced with corresponding annihila- 
tion operators and complex conjugates — with Hermitian 
conjugates. 

We represent the optical field inside the cavity as a 
sum of two light waves: the one A + (x,i) running in the 
positive direction of the x-axis and another one A-(x, t) 
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running in the opposite direction. We write them down 
as sums of three parts. The first part is the "strong" 
wave with amplitude A±q. The second part is the ex- 
plicit formula © for direct coupling 'GW+EMW. The 
third part is the "weak" addition a±(x, t) — the unknown 
function of its arguments representing both optical noises 
and phase shift acquired while circulating inside the cav- 
ity. 

A±(x,t) =A ±0 e- i ^ otTkox) +A ±0 g±(x,t)e- l ^ otTko ^ 
+ a±(x, t) e -<(«o*T*o») + c . c-j ( lla ) 



a±(x,t) = 



o±(n + uj )e~ 



iO(t=Ff) 



dn 

2?r' 



(lib) 



We assume interaction between a gravitational wave h 
and "weak" field a± to be a negligible effect since (i) it is 
most likely that coupling of GWs to optical noises does 
not contain any useful information and (ii) coupling of 
GWs to phase shift (variation of optical field) acquired 
due to GW itself leads to the negligible h 2 terms. 

Output wave A out (x,t) includes the wave reflected 
from the input mirror and the wave flowing out of the 
cavity. The latter one carries information about grav- 
itational wave (signal wave) along with optical noises. 
Again we represent the wave as a sum of "strong" and 
"weak" components: 

A out (x,t) =A ouW e-^° t+k ^ + a out (M)e- i(< "° t+fcox) 



+ C.C., 

r+ao 

dout(x, t) = / a ou t 



(f! + w )e- 4O ('+t) 



2w' 



Our problem is to find connection be- 
tween A ±Q , A ouW , a±(x,t), a out (x,t) and 
A- m o, di n (x,t), g±(x,t). This connection is derived 
from the boundary conditions. The first boundary con- 
dition is the requirement of vector potential continuity 
along the input mirror: 

A out (0, t) = RA in (0, t) + TA_(0, t), (12) 
A+(0, t) = TA ia {0, t) - JL4_(Q, t). (13) 

The second condition is the requirement of electric field 
to be equal to zero on the surface of end- mirror in its local 
reference frame [Hj]: E' z (0, t — r) = (we neglect X 2 /c 2 - 
terms), where notation r = L/c is introduced. Using 



well-known formulas for Lorentz transformations of elec- 
tric and magnetic fields from movable frame of reference 
to the laboratory frame {33] we rewrite this condition in 
the following form: 

E Z (L + X, t - t) + (X/c)H y (L + X,t-r) = 0, (14) 



or, substituting here definitions E z 



[l/c)(dA/dt) and 



-dA/dx, we come to the following equation: 



dA(L + X,t-r) -dA{L + X,t-T) 

h X 



= 0, 



dt dx 
which means that directional derivative of vector- 
potential equals to zero along the mirror trajectory. Thus 
it is sufficient to demand A(L + X, t — r) = or 



A+{L + X, t - t) + A_(L + X, t - t) = 0. 



(15) 



Here X = X(t — r). We would like to stress that this 
boundary condition automatically allows us to take into 
account the terms proportional to X/c which lead to the 
appearance of the force of radiative friction (see Sec. IVI|) . 

Strictly speaking, we cannot use Lorentz transforma- 
tions to change the reference frame in our case since 
space-time is curved. However, one can verify that cor- 
rections to formulas of Lorentz transformations due to 
gravitational wave are negligible. 

We assume further that mirror's displacement X{t) ~ 
Lh <C L. Splitting the set of equations (fT2l [T3l [T5"|) into 
the zeroth and the first order sets we obtain their solu- 
tions. It should be noted here that we have not taken 
into account the influence of GW on time in boundary 
conditions. For the zeroth order this effect vanishes since 
dependence on time is only included in the phase factor 
e -u*j t which i s a common multiplier for both the left and 
the right sides of equations. For the first order this effect 
is negligible since it will lead to terms of the second order 
in h. Therefore, the zeroth order solution is: 



.4 



T 



1 - Re 2iuJ ° T 



A 



inO: 



.4, 



outO 



1 - Re 2iuJ ° T 
R - e 2iu,0T 
1 - Re 2lu> ° T 



(16a) 
(16b) 



A 



inO, 



and the first order solution in spectral domain is: 



J 



a + (tt + CJ ) = Oin(f} + ^o) 



T 



Re 2lulor A +0 



gw+emw (n) 



(17a) 



a_(fi + lu ) = — ai n (0 + lu ) 



A-o 



n nvv 

(fi) 



(17b) 
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<W(^ + Wo) = «in(^ + W ) 



I? 



,2i(O+w )t 



1 _ Jl e 2i(n+uj )T 



TA_ 



1 — i2e 2< ( n +"'o)T 
I 



j* mir (n) + 5* 



gw+emw 



(Q) 



(17c) 



Here we introduced the convenient notations for phase 
shifts after single reflection: 



5V mil (n) = 2k X(Q)e 



iQ,7 



(18) 



is the phase shift acquired by the light wave when reflect- 
ing from the movable mirror. Note, that in time domain 
it comes with delay t — r since it takes light wave time r 
to travel from the end-mirror to the input mirror; 



8f 



(ft) =i g- (L, ft + lu ) — g + (L, Q + u) ) 



gw+cmw ^ 

is the phase shift acquired by the light wave from its 
direct interaction with gravitational wave. Below we cal- 
culate it explicitly and show that in long-wave approxi- 
mation this phase shift vanishes. 



V. PHASE SHIFT OF LIGHT WAVE IN THE 
CASE OF SINGLE REFLECTION 

Let us compare our results with the ones already de- 
scribed in literature for the case of single reflection. To do 
this we ignore all resonant features of our system. First, 
we shall rewrite our results in spectral domain explicitly. 
For direct coupling of gravitational wave to light wave 
(Un]), taking into account w gw <C ujq, we get: 



gw+cmw 



(Q) PS -k Lh(Q) 1 - 



sin f2r 



(20) 



In time domain we consider the case of a monochromatic 
GW with spectrum h(Q) — 2Trh S(fl — cj gw ) + 2irh 5(fl + 
). After simple integration of formula (j2"0)) we obtain: 



gw+emw 



(t) = — 2k Lh cosu gw (t — t) 



ho 



sinWg W < — sinw gw (i — 2r) 



It is straightforward to verify that in long-wave approxi- 
mation L/A gw <C 1 this formula reduces to 



6V 



gw+cmw 



(*) 



(2tt) 2 , T1 ( L . 
- koLriQ y- — ) cosw gw t 



+ o 



k Lh ( - — 

Ag W 



Let us now turn to the full phase shift 5 'I'll — 



gw+emw ■ 



Until we find mirror law of motion we do not 



know how displacement X(fl) depends on gravitational 
wave h(£l). However, for single reflection (the case when 
we can ignore the force of light pressure) this solution has 



been known long before: X(Q) = ^Lh(fl). Therefore, 
the full signal is 



S^ LL (Q) = k Q Lh(Q)^^e^. 

ilT 

In time domain for monochromatic GW we get: 



(21) 



^LL(t) - ho 



to 



sinwgwt — sinu; gw (i — 2t) 



This result coincides with the one obtained in the frame- 
work of the TT gauge in Ref. [5|: 



TT 




h{t')dt'~ \ h(t')dt' 



sin uj gw t — sino>g W (£ — 2t) 



LL- 



This coincidence confirms the results of Rcf. [2l|: for- 
mulas for the net phase shift obtained in the TT and the 
LL gauges are equivalent if only in LL gauge one takes 
into account the effects of localized and distributed grav- 
itational redshifts (along with displacement of mirror). 
Remind that distributed redshift is the direct interac- 
tion of gravitational wave with light wave. Note that 
we have not considered the effect of localized redshift 
simply because our clocks are located at x — where 
qori^O, t) = — 1. For generalization one may turn to Ref. 

Again it is straightforward to verify that in long-wave 
approximation the full signal reduces to 

L 



S^LL(t) pa 2koLho cosw gw i + O I koLho 



A 



thus leading to the rising of the term containing 
(L/A g w) ho in the full phase shift. This was not men- 
tioned in Ref. [j| since the delay t — r of the reflected 
light wave was not taken into account in J^mh- there: 
£*mir(0 = 2k X(t - t) but not 2k a X(t). 



VI. MIRROR EQUATION OF MOTION 

In the absence of non-gravitational forces test mass 
moves along the geodesic (with respect to the input mir- 
ror of the cavity where the origin of coordinate frame 
is set) as dictated by geodesic equation [f| [33j]. It is 
straightforward to verify that Christoffel coefficients are 
Too = —xh/2c 2 , = = (for the motion along 
the x-axis) in metric @, so in non-relativistic limit we 
obtain the following equation of motion for a free mass: 



fx 

IF 



i 



xh{t) = 0. 



(22) 
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We stress here that this equation is strict for any sepa- 
ration between the test masses. Its solution is trivial if 
obtained with the method of successive approximations. 
By assuming x{t) = X^\t) + X^{t) with = L wc 
get =X{t) = \Lh(t). 

Generalization of Eq. (f2"2"]l for the case of non-geodesic 
motion is the following: 

d 2 x 1 ; - , . 1 ^ . . . 

-T7T _ o Xh ^> = _ Fx ( X > X > *)' 

dt z 2 Tji 

where F is the force of non-gravitational nature, the force 
of radiation pressure in our case, and m is the test mass. 

In general, the force of light pressure should be cal- 
culated in curved space-time by means of General Rel- 
ativity. However in Appendix [B] we show that this cal- 
culation can be performed in the framework of classi- 
cal flat space-time electrodynamics. We should explain 
this statement in detail. The rule for calculation of the 
light pressure force in flat space-time for any given elec- 
tromagnetic wavefunction A(x,t) is known [331 ] - Obvi- 
ously, there exists a general rule for calculation of light 
pressure on curved space-time background. Our state- 
ment is that this general "curved space-time rule" dif- 
fers from "flat space-time rule" only in negligible cor- 
rections. Therefore, in our case we can apply the "flat 
space time rule" to A(x, t) which includes gravitational 
wavefunction h (via direct coupling of GW to EMW) 
without introducing any significant error. Namely, in 
Appendix [5] we estimate the resulting inaccuracy as 
dX^l = F^l/mQ 2 « {£ FP /mc 2 )Lh{t) « lQ~ 17 Lh(t) 
(here £ fp is the full energy of light wave inside the cav- 
ity) compared to the first order solution X^ — | Lh(t) 
of Eq. (|2"2"|) , assuming that £fp = 20 J and m = 40 kg 
(typical parameters for Advanced LIGO project). 

In classical electrodynamics [33| spatial components of 
the radiation pressure force are F % — J T lk dSk , where 
T lk are the components of the energy-stress tensor of 
electromagnetic field and dSk is the oriented area el- 
ement. Since the light beam with cross-section S is 
parallel to the x-axis and normal to the mirror sur- 
face then the x-component of the light pressure force 
is F x = T n S = T 1X S, where T n = (E 2 + H 2 )/8n, 
E z = -(l/c)(dA/dt) and H y = -dA/dx. Therefore, 
we obtain the following first order equation of motion: 



d 2 X 1 •• , . 
-W-~2 Lh ® 



= F^(X,X,t), 



(23) 



where 



F^(X,X,t) = 



S_ 

8tt 



OA 
~dt 



OA 
dx 



(i) 

x=L+X(t) 



Index (1) on the right side of equation means that only 
the first order in h and X terms should be considered. 
Note that the zeroth order force of light pressure is con- 
stant, therefore, we should redefine the length L in all 



formulas, assuming that in real experiment constant force 
is compensated somehow. The first order force is calcu- 
lated explicitly in spectral domain in Appendix [Cl (below 

we use notation F x instead of Fx )■ 



F x {n) = F Ruct (n) + F pm (n) + f ( 



gw+cmw 



where 

-Ffluct(^) 

Sk 



T 2 A* 



(lo + fi)e 4 



f 



a in (w + ft) 

I _ R e 2i^ T I _ Re -2i(uJo-n)7 a m( W - ft) 



Re~ 2iuJ or i 
T 2 A in0 



(lu - ny nr 



is the fluctuation part of the force (see Appendix IC 2p , 
which leads to the quantum radiation pressure noise 
(back-action) after the quantization of the light field; 



Fpm(ft) 



/C(ft)X(fi) 
K(n)X(Q) 



2iar(n)x(n), 



is the ponderomotive force, where (see Appendix IC 3P : 



4k SW FP Re MtT sin2u; T 
1 - 2Re 2inT cos 2uj t + R 2 e 



UQ.T ' 



is the coefficient of optical rigidity and 



r(ft) 



SWi 



FP 



l-JFe 



2MQ.T 



1 - 2Re 2lUr cos2uj T + R 2 e 



-UV.T 



(24) 



(25) 



is the coefficient of radiative friction. Here Wp-p = 
kl(A +0 A* +0 + A- A*_ )/2tt is the energy density of light 
wave circulatin g in side the cavity. Formula (|24[) was ob- 
tained in Refs. [IlElEiJ] and formula ((25]) in Ref. [jj], 
where phenomenon of radiative friction was discussed in 
detail. 

We split up JC into K and 2iQT using the following rea- 
soning. The obtained mechanical impedance Z(f2) of the 
mirror is Z(Q) = mfl 2 +2inr(Q) — K(fl) which obviously 
looks similar to the impedance of linear harmonic oscil- 
lator. Therefore, /C(ft) includes the term 2«f2r(f2) de- 
scribing friction, and the term K(Vt) describing rigidity, 
where 2iQT is of the order of (Q,/ujq)K . Radiative friction 
r(f2) appears due to the boundary condition (fT4"]) which 
takes into account the terms proportional to X/c. This 
relativistic correction was not considered in Ref. . 

Both the terms have real and imaginary parts. The 
latter means that both effects are introduced with de- 
lay. Note that the coefficient of "true" friction is only 
the real part of formula (f2"5]l . "True" friction is always 
masked with the imaginary part of formula (|24p since 
| if | ^> |2ftr|, except for the case of optical resonance (see 
later) when optical rigidity vanishes. In non-resonant 
regime there exist two possibilities (in long-wave approx- 
imation), depending on the sign of sin 2ujqt: either (i) 



8 



m[K{n)] > (positive rigidity) and 3[-KT(f2)] < (neg- 
ative damping), or (ii) 9l[if(il)] < (negative rigidity) 
and 3[if(fi)] > (positive damping). These inequalities 
were first obtained in Refs. 0, @j- 

The concept of ponderomotive rigidity and the role 
it plays in standard quantum limited (SQL) detectors 
has been widely discussed in literature before (for long- 
wave approximation only) . For the first time it appeared 
in Refs. 0, [|| where microwave resonators were consid- 
ered. Optical ponderomotive rigidity and its application 
to large-scale interferometric GW detectors were studied 

in Refs. @, ES EL El- Papers 0, EI EE El E3 

dealt 

with the optical springs in signal-recycling topologies of 
LIGO detectors; 



*gw+cmw(^) — 2 ^(^) ( 1 



sin fb 



Lh(Cl), (26) 



is the 'GW+EMW correction to the force of ra- 
diation pressure (see Appendix IC 4j) . arising when 
the cavity length becomes comparable with gravita- 
tional wavelength. Note, that it is proportional to 
<5^E'g W+emw (f2) much as ponderomotive force is propor- 
tional to <5\& m i r (CI) . Obviously, in long- wave approxima- 
tion L < A gw or f2r < 1 this 'GW+EMW correction 
vanishes (along with the corresponding phase shift). Sim- 
ilar formula (without T in K) was obtained in Ref. [ill ]. 



VII. MIRROR'S LAW OF MOTION AND 
DETECTOR'S RESPONSE SIGNAL 



Now we solve Eq. (|23f in spectral domain: 

1 mfl 2 Lh(fl) F gw+cmw (n) + F fluct (fi) 



X(Q) 



2 mf2 2 - JC(O) 



in 2 - jc(o) 



(27) 



In this paper we are mostly interested in dynamical part 
of the problem, therefore, we shall ignore fluctuations 
(terms containing a m and Ffl uct ) further. Substituting 
(l26l) into J27D we obtain: 



X(Q) 



Lh(Cl) 



ifl 2 - K(Cl) 



mfT-/C(fi) 1 



sin f2r 



(28) 

Strictly speaking, in order to keep all the X/c- 
proportional terms we should expand g± in formula 
(fT9|l up to the first order of f2/u;o and then replace 
(1 — sine fir) with (1 — fi/wo)(l — sine fir) in formulas 
([26]) and (|2SJ). However, for S1/2tt - 10 2 10 5 Hz we 
have fi/wo ~ 10~ 12 -j- 10~ 9 and this correction can be 
omitted in both the cases K (f2) ^ (term K x fl/u>o is 
negligible) and K(Q) = (term 2ttT x Q/u oc X 2 /c 2 
is negligible). In addition, the term in square brackets 
in formula (|28p does not contain corrections of the order 
of (£fp / mc 2 )m£l 2 « 10 _17 mri 2 since we have omitted 
them when calculating the radiation pressure force. 



Substituting formulas (J20J) and (28J into (|17c[) we fi- 
nally obtain the explicit formula for the optical field of 
the signal wave: 



Oout(^ + Wo) 



(1 - Re 2iuJ " T ) (l - Re 2l ( n +^o)r) 
mfl 2 , _ . sin fir 



in 2 - JC(fi) 



ikoLh(Q) 



fir 



(29) 



Similarly to the previous formula one can replace sine fir 
with [l — (1 — Q,/l>o)(1 — sine fir)] for strictness, but for 
the reasons discussed above this correction is negligible. 

Note that the obtained response signal is the product of 
three terms: optical resonant multiplier, mechanical res- 
onant multiplier and the term (|21[) describing the single 
reflection response. The latter is proportional to sinc- 
function and is sometimes referred to as the "electric" 
component of the detector full (single reflection) angu- 
lar pattern [3||. Recently it was demonstrated that the 
angular pattern also contains the so-called "magnetic" 
component [Ii| (due to non-optimal GW polarization 
and detector orientation in general case) with its am- 
plitude proportional to Or in low frequencies but rising 
with increasing frequency [361 ]. Therefore, it may be in- 
teresting to perform a (straightforward) generalization 
of our results for the case of generic GW polarization 
and direction of propagation. Probably, it will add up 
to the multiplication of the known single reflection angu- 
lar pattern by the same optical and mechanical resonant 
multipliers as in formula (|29[) but, strictly speaking, this 
problem requires a detailed analysis. 

In many cases the phenomenon of radiative friction 
is negligible in comparison with optical rigidity since 
|20r| < K. Then one should replace fC(Q) -> K(Q) 
in formulas j28j and |29|) . 

One of the most interesting effects which may arise 
when detecting the high-frequency GW signal is para- 
metric transitions between optical modes of the cavity 
if Q = uj gw « 7iajf sr , where n G N and u% r = nc/L is 
the cavity free spectral range (see Fig. [3]). In some sense 
this effect looks similar to the phenomenon of parametric 
oscillatory instability (POI) described in Refs. [37l. [38|. 
The difference is in the following: POI appears due to 
the coupling of optical modes to acoustic (mechanical) 
modes of the movable mirror while the effect illustrated 
in Fig. [3] appears due to the coupling of optical modes to 
gravitational wave itself. This effect can be used for de- 
tection of high-frequency signals and requires additional 
analysis. 



VIII. SPECIAL CASES 

It is useful to examine the obtained results in two spe- 
cial cases — optical resonance and small detuning from it 
— and compare our results with the ones already present 
in literature. We shall use notations for the energy of 
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FIG. 3: Optical resonances of Fabry-Perot cavity: uj c is one 
of the eigenfrequencies, Wf sf = nc/L is the cavity free spectral 
range, 7 = (1 — R)/2r is cavity half-bandwidth. Red dashed 
arrows denote some possible transitions between neighboring 
resonant curves under the influence of GW with frequency 
equal or approximately equal to FSR: oj gw ~ Wf Br . Blue dash- 
dotted line denote one of the possible transitions between non- 
neighboring curves if u) gw ~ 2uf sr . 



light wave inside the cavity £pp = W F pSL, and cav- 
ity half-bandwidth 7 = (1 — R)/2t. In these notations 
T 2 = 1 - R 2 « 4 7 r and £ FP = 2P in / 7 , where P in is the 



input power. 



Resonant features of the cavity only appear in frequency- 
dependent amplification of the light wave amplitude. 

In long-wave approximation these formulas, denoted 
with 'lw', reduce to the well-known ones: 



pros 
1 lw 



1 j ~ l 2 j -in' 



The coefficient of "true" friction is the real part of r[^ s : 



«[r;-; s (^)] = -ji 



£ FP 7 



32P ir 



1 



L 2 7 2 + ft 2 c 2 T 4 + (4fir) 2, 
Last formula with P m coincides with the one in Rcf. It 



1 



x^(n) = -Lh(n) 



1 



in 2 + 2inr™ s (n) 2 



Lh(n), 



a out, lw(^ + ^0) 



1 



-A; 



(7 - £0)t mn 2 + 2ifl rf° s (fi) 

1 



ikoLh(fl) 



inO" 



ikr,Lh(Q). 



(7 — iVl)t 

B. Small detuning from optical resonance 



A. Optical resonance 

Let us assume that frequency of the pump wave loq co- 
incides with one of the eigenfrequencies of a Fabry-Perot 
cavity: uioT = nn, n E N. We shall denote all quan- 
tities with index 'res' in this section. Then, simplifying 
formulas flU}, (|2"5|). (|2"5]l and (QU) wc obtain 



K les (fl) = 0, r res (ft) 



SW F p 1 + Pe 



2if2T 



l-i?e 



HVlT ■ 



thus in the case of resonance rigidity vanishes and radia- 
tive friction is not masked by it; 



X res (Q) = ~Lh(Q) 



1 - 



2£nr res (f2) sin Or 



mn 2 + 2ifir res (fi) fir 



1 - i?e 2inr mfi 2 + 2inr res (rt) 



x ik Lh(fl) 



sin fir 
fir 



However, radiative friction is still a tiny effect even in the 
case of optical resonance and is omitted usually: 



X re8 (0) a i L/i(0), 



l-i?e 



2if2i 



■ ik()Lh(fl) 



sin fir 



Therefore, in the case of optical resonance light field in- 
side the cavity does not act dynamically on test mass 
if we neglect radiative friction, i.e. test mass is inertial. 



We shall now turn to the case when frequency of the 
pump wave is slightly shifted from one of the eigenfre- 
quencies: 

uj a = \-6, 6t<1. 

T 

Again, simplifying formulas (|24p , [)28p and (12T))) we obtain 
(all quantities in this section are marked with 'det'): 



K de \rt) 



8k SW FP Re 2lUT ST 



(1 - Re 2inT ) 2 + 4i?(<5r) 2 e 2iOT ' 



(30) 



The graphs of 9\[K dct (n)] ("true" rigidity) and 
3[K det (Q)] (damping) are plotted in Fig. @. We do not 
consider radiative friction here since ^> |20r(fi)|; 



x dct (n) =-Lh(n) 



K de \n) 



O ^(n + wo) = -^ 



in() 



sin Qt 
' mfl 2 - K dct (Q) VLt 
27 1 + 2iSr 



j-iS 1-R(1 + 2iST)e 2inT 
"'" l k Lh(tt) Sh,i}7 



in 2 - K det (n) 



Or 



(31) 



Therefore, optical rigidity appears only when pumping 
frequency is detuned from resonance. 

In long-wave approximation ('lw') we obtain: 



FP 



L 2 8 2 + { 1 - iQ) 2 



(32) 
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FIG. 4: "True" rigidity (upper figure) and damping (lower 
figure) normalized to their maximum values. For visualization 
we assumed R = 0.6, 5/2n — 100 Hz. 



coinciding with the corresponding formula in Ref. 
One can easily check that the ratio 9\[K det (fl)]/(mfl 2 ) 
can be made of the order of unity for the values of param- 
eters planned for Advanced LIGO (L = 4 km, £pp < 20 
J, to = 40 kg, lu /2tt = 3 x 10 14 Hz, 7 /27r « 1 Hz, 
w 5/2tt = 100 Hz); 



Xdet 
lw 



-i/i(0), 



2 out, lw(^ +^o) 



mft 2 - ^ d w ot (0) 2 

7/r 1 + 2i<5r 



(33) 



A 



in() 



-7 — i<5 7 — i(<5 + 0) 
mfl 2 



mfl 2 



Z^dct 



(fi) 



ik Lh(Cl). (34) 



Formulas (|32H34j) are applicable to consideration of the 
signal-recycling topology since it is equivalent (with some 
modifications of parameters) to a single detuned Fabry- 
Perot cavity [ll|. Therefore, all the major results of Refs. 



[131 . Il4 116L |17| which studied in detail optical springs in 
SR topologies hold true for a single cavity. 



C. Response function of the detuned cavity in the 
vicinity of FSR 

To analyze the operation of GW detector near the FSR 
we turn to formulas (|3"0"|) and (|3"Tj) . Introducing notation 
A = fl — Wf sr , |A| <C Wf sr we rewrite these formulas in the 
following way: 



K det (A) = 



fp 



L 2 5 2 + (7 - iAy 



in the vicin- 
u>{ sx ) in Fig. 



We plotted <R[K det (A)] and ?[K det (A) 
ity of A = (or, equivalently, near Q 
[H Note that while decreasing the cavity half-bandwidth 
from 7 ~ 5 to 7 <C 6 the shape of $H[i4r dct (A)] changes 
from "singlet" to "doublet". Obviously, this takes place 
in the vicinity of each f2 = najf sr , n £ Z; 



2 out ( W fs: 



.4 



inO ' 



A + LU ) 

7/r 



1 + 2iSr 



7 



i5 7 

J fsr 



i(S + A) 
A 



^fsr 



K det (A) uj ft 



ikoLh(uj{ sr + A). 



The ratio m[K det (A)] / (mujfj can be estimated as « 
10~ 4 (see Fig. © for L = 4 km, w /2tt = 3 x 10 14 Hz, 
r5/27r = 100 Hz, 7 /27T = 1 Hz, £ FP = 20 J, m = 40 
kg and Wf sr /27r = 37.5 kHz. Thus the implementation 
of optical rigidity for the detection of signals near the 
FSR will require the increase of circulating optical power, 
the decrease of cavity half-bandwidth (see Fig. [S]) and 
the mirror mass. With these requirements met it will 
be possible to perform the resonant detection of high- 
frequency gravitational waves in the vicinity of 37.5 kHz. 
Note that while mu)f sx — 9\[K dct (A)] can be made vanish 
for certain frequencies in principle, the resonant gain will 
still be limited by mw| r /3[^ dct (A)] . 



IX. CONCLUSION 

In this paper we have analyzed the interaction of '+'- 
polarized plane gravitational waves with a Fabry-Perot 
cavity. Calculations have been performed in its input 
mirror's local Lorentz frame which is best suited for 
consideration of non-gravitational forces acting on test 
masses, such as the force of radiation pressure. We have 
considered the general case of ratio L / A gw (L is the cavity 
length) for the purpose of high-frequency GW signals de- 
tection when long-wave approximation is violated. Such 
high-frequency gravitational waves can be emitted by the 
sources predicted in modern astrophysical and cosmolog- 
ical theories, for example, relic GW background. 
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FIG. 5: Normalized rigidity (upper figure) and damping 
(lower figure) as functions of detuning from the FSR. For 
visualization we assumed 5/2ir — 100 Hz. 



FIG. 6: Rigidity (upper figure) and damping (lower figure) 
near the FSR compared to mu^. For visualization we as- 
sumed L = 4 km, lo /2-k = 3 x 10 14 Hz, 5/2-k = 100 Hz, 
£fp = 20 J, m = 40 kg and w fsr /27r = 37.5 kHz. 



We began from the solution of wave equation for run- 
ning light wave which is coupled to a weak plane grav- 
itational wave and found that in long-wave approxima- 
tion this interaction leads to a small correction of the 
order of (rr/A gw ) 2 h, where x is the distance traveled by 
light, in agreement with Ref. However if x ~ A gw 
direct interaction between gravitational wave and light 
wave becomes significant and is described by oscillating 
function. We applied this result to consideration of the 
boundary problem for the light wave in a Fabry-Perot 
cavity with the movable mirror. It was shown that the 
response signal of such GW detector includes two terms: 
(i) the phase shift of light due to reflection from the mov- 
able mirror and (ii) the phase shift due to direct coupling 
of gravitational wave to light wave. In our analysis we 
assumed that the end-mirror of the cavity moves under 
the influence of (i) gravitational wave and (ii) the force 
of light pressure, neglecting all other forces. As a special 



case, by neglecting the resonant features of the cavity, 
we examined the situation of single reflection when radi- 
ation pressure is negligible. In this case test mass under- 
goes geodesic motion only and, therefore, is inertial. TT 
gauge can be applied then to the analysis of GW detec- 
tor. This was performed in Ref. [5] and we found that 
our LL-based results coincide with the TT-based results. 
This also confirmed the results of Ref. [2l[ : in order to 
get the agreement between the LL- and the TT-results 
direct interaction of the gravitational wave with the light 
wave must be taken into account in the LL gauge along 
with the displacement of test masses. 

Next we considered the mirror equation of non- 
geodesic motion in the field of gravitational wave and 
the force of light pressure. Traditionally in literature 
this is done only in long-wave approximation. We pro- 
posed rather general method of calculation of radiation 
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pressure when condition L <C A gw is violated and showed been assigned LIGO doc. number LIGO-P060055-01-R. 
that this calculation can be performed in flat space-time 
without introducing a significant error. We obtained the 
formulas for coefficients of ponderomotive optical rigidity 
and radiative friction along with the general relativistic 
correction describing the influence of direct 'GW+EMW 
interaction on the force of radiation pressure. This cor- 
rection must be taken into account in order to get correct 
formula for the response signal when long-wave approxi- 
mation is violated. Though formulas for optical rigidity, 
radiative friction and 'GW+EMW correction separately 
have been presented in literature before [ill [l2l.Tl8l]. no- 
body have obtained them simultaneously using one uni- 
fied method, as we have done in this paper. The phe- 
nomenon of optical rigidity plays significant role in op- 
eration of proposed schemes for advanced GW detectors 
and our results may be of potential interest in consid- 
eration of the regimes of high-frequency GW signals de- 
tection. One may also apply our results to calculate the 
modification of the detector angular pattern, including 
its magnetic component [35,36], due to radiation pres- 
sure effects. 

The possibility of parametric transitions between dif- 
ferent resonant curves of the cavity under the influence of 
gravitational waves with frequencies u> gw w nu>[ SI was dis- 
cussed. This effect may find implementation in detection 
of the high-frequency GW signals using optical modes 
with frequencies ujo ± nuj{ SI as signal waves instead of the 
wo-input mode. 

In the end we considered two special cases of optical 
resonance and small detuning from it. Excellent agree- 
ment with results of [HI, [HI was found. We also esti- 
mated the deposit of optical rigidity into the response 
signal near the FSR and showed that in order to per- 
form the resonant detection of GWs near that frequency 
one should increase the circulating optical power and de- 
crease cavity half-bandwidth and the end-mirror mass. 
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APPENDIX A: SOLUTION OF THE FIRST 
ORDER EQUATION FOR RUNNING EM WAVE 

In this Appendix we obtain the solution of the Cauchy 
problem (JSHZJ)- 

It is convenient to go to spectral domain. Applying the 
theorem of convolution to the right side of Eq. ((5|), we 
get the following first-order equation in spectral domain: 



d 2 A^ fi 2 
9a; 2 c 2 



-k 2 



A il) 

A+0 e^ h(n coo) 



A%e - lk0 A^±pl h{n + LOQ) 



+ A- Q e- ikoX 
+ A 



(Q - w ) 2 



h(Q — loq) 



elkoX (» + w ) 2 h(n + UJo) 



where 



AW(x,n) = A^(x,n) + A^\x,n), 
A i + ) (o,n) = aL x) (o,q) = o. 



(Ala) 

(Alb) 
(Ale) 



are equivalent to formulas © and ([7]) correspondingly. 
For convenience we introduce the following notations: 



A^(x,Cl) = A 



— kr,A 



B 2 
B 3 



1 



= th 2 A* 
— 2 ^o^+o 

= 5 k l A -* 



1 



k 2 A* 





k = 


(fi 


-wo) 2 




c 2 


(0 


+ w ) 2 




c 2 


(CI 






c 2 




+ ^o) 2 



n 

c 

h(Cl — too) 
h(Q + ujo) 
h(Q — luq) 
h(Q + ujq). 



2 - • c 
Eq. (|Ala[) rewritten in these notations is 

cPA 
dx 2 



-k z A = (Bx+B^e^ + (B 2 +B 3 )x 2 e 



-ikox 



(A2) 



Note that spectral frequency £1 (or k — fl/c) plays the 
role of parameter here. General solution of homogeneous 
equation d 2 A/dx 2 + k 2 A — is well-known: A(x) = 
C ie ikx J r c 2 e~ lkx . We shall obtain the particular solution 
of heterogeneous Eq. (|A2j) using the method of variation 
of constants (see any ODE handbook), assuming 



A{x) = d{x)e lkx + C 2 (x)e 



-ikx 



(A3) 



Following the general rule of this method we must solve 
the following set of equations: 



dCi e ikx 
dx 



dQ 
dx 



2 ^—ikx 



=0. 



ik dCi e , kx _ lk <^l e -«x =(B + Bi )x 2 e ik ° x 
dx dx 

+ (B 2 +B 3 )x 2 e- lkox . 
It is straightforward to verify that solution is 



Cl (x)=^±^e^ k "-^ x 



2x 



j{k - k) (fc - k) 2 i{k - k) 3 
B 2 + B : 



2ik 
X ( — 



3 e -i{k a + k)x 



2.r 



i{k Q + k) (fco + fc) 2 i(k + k) 3 



+C-_ 



in- 



and 



C 2 (x) =- Bl + B \ ^o+k)x 



2.x 



i{k + k) (fc + k) 2 i(k a + k) 3 



B2 + B 3 e ^ ko -k)x 
2ik 



2x 



i{k -k) (fco-fc) 2 i{k -k) 3 



-C: 



20- 



Substituting these solutions for C\(x) and C 2 (x) into 
formula (|A3|) and keeping in mind our notations and for- 
mulas (|9aH9cl [Alb|l we obtain the general solution for 
"positive" and "negative" waves: 



A^'{x,fl) = C ±0 e 

+ l -A ±0 h{Vl-LO )e ±lkoX 

" x 2 ^(VL) =f ixrj(Q) -((fl) 

+~A* ±0 h(n + u>o)e* ik ° x 

(-fi) ± ixri(-Sl) - C(-O) 

(A4) 

Here we redefined for convenience constants C10 and C20 
as C + o and C_o to describe separately the additions to 
"positive" and "negative" waves. Constants C±o are ob- 
tained from initial conditions (jAlcP : 

c± = \ A ±0 h(n - wo)c(n) + \ A* ±Q h(n + wo)C(-n)- 



Substituting them into formula (|A4j) and performing in- 
verse Fourier transformation we finally obtain the first 
order solution in time domain (El). 
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APPENDIX B: GENERAL RELATIVISTIC 
CORRECTIONS TO THE FORCE OF LIGHT 
PRESSURE 

In this Appendix we shall prove the statement made 
in Sec. I VII the 4-force of light pressure Tg m can be 
calculated by means of classical fiat space-time electro- 
dynamics without introducing a significant error. 

Whatever the general relativistic rule for calculation of 
the light pressure force in any case it will be proportional 
to: (i) the components of energy-stress tensor T a ° of the 
light field inside the cavity averaged over period of optical 
oscillations 2ir/ajQ, (ii) cross-section of the light beam S 
and (iii) components of metric tensor (J2J: 

ex g Xv T a PS. 

The energy-stress tensor of electromagnetic field in 
curved space time is [HI: 



1 

47T 



F F aa a P 



, r af3 r yfiv 



where F a p is the tensor of electromagnetic field A 1 inside 
the cavity: 



F a p =V a Ap - VpA a = gp^VaA 1 - g ai VpA J 
' dA 1 



dx a 



r 7 A v 



( dA'< 

9aj [d^ 



r 7 a v 



Here V Q is the covariant derivative with respect to the 
connection of space-time r^ A : 



r' 



1 „«/ dg av dg D 



dg„> 



dx v dx° 



We shall now examine the orders of the values included 
in F a p. To do this we shall omit all numerical constants 
of the order of unity and all terms of the second and 
greater orders in h. Besides we shall consider only the 
case of monochromatic GW with frequency w gw since 
generalization is trivial. Remind that all functions are 
evaluated at point x = L. Returning to Sec. U we con- 
clude that 



1 



I? 



1 v\ 



T 2 

72 



c 



L 



Remind that (see formula (|1 la|) electromagnetic field 
A 1 {x,t) inside the cavity is a sum of "positive" wave 
A + (x,t) and "negative" wave A-(x,t). These waves 
include "big" components with amplitudes A±q and 
"small" components proportional to A±oh (we are not 
interested here in the term proportional to X). It is 
necessary to mention that these "small" terms are in- 
cluded in the wavefunctions of electromagnetic field and, 
in general, the gravitational nature of these small cor- 
rections is of no importance to us here. Therefore, we 
shall not consider them as general relativistic corrections 
further, since we are interested only in GR corrections to 



the methods of flat space-time electrodynamics but not 
the corrections included in A 1 {x,t) itself. We estimate 
the order of the derivative of EMW field as 



dA 

dx 



k Ao + k A H{L,t), 



where H(L,t) ~ h is the mentioned small correction to 
the EMW field inside the cavity, containing GW wave- 
function. Thus we can estimate the order of F a g in the 
following way: 



F aB ~ k A Q ( 1 + H(L,t) + % h + ^ ^ h + A 1 



0J 



UJ C 



Here Aq has the order of magnitude of the electromag- 
netic wave's amplitude inside the cavity. The light pres- 
sure 4-force T£ m is proportional to F a p, "squared" (re- 
mind that we neglect all numerical coefficients), averaged 
over period of optical oscillations 2t:/luq (terms propor- 
tional to e 2lw °* disappear) and multiplied by S: 

J^ m ~W FP S 

n + H(L,t) + ^h + ^^h + — h 



UJ Q C 



where W^fp ~ ^qAq 1S proportional to the energy density 
of light wave inside the cavity. 

Let us now return to the mirror's equation of motion: 



d 2 x^ 



^ dx" dx x 
vX ~ds~ds 



1 



We rewrite this equation in non-relativistic limit for x- 
component, expanding the force of light pressure into the 
zeroth and the first order in h terms: 



fx 



where 



f(°) ~ W FP S, 



F, 



GR ) ' 



is the constant force of light pressure; 

F$ ~ W FP SH(L,t), 

is the part of the force proportional to the phase shift of 
light due to its direct interaction with GW; 



F 



Lr 

c 2 



£i~W FP S[^h+^^h-r—h 



L 2 v 

6 



L 

U C 



is the GR correction to the methods of flat space-time 
electrodynamics. Our problem is to prove that it can be 
omitted when calculating the force of light pressure. We 
use the method of successive approximations to solve the 
equation of motion. The zeroth order force is constant, 



15 



therefore we should redefine the length L in all formu- 
las, assuming that in real experiment constant force is 
compensated somehow. The first order equation is then 



d 2 X 1 



dt 2 



Lh(t) 



■F 



(i) 



FP 



1 

m ' 



To estimate the orders of all the terms included in this 
equation we turn to the spectral domain: 

-si 2 x{si) + ^Lii 2 h(n) = —F$(n) + — f£1(Q). 



Its solution is: 

x(si) ~ Lh(n) 



W FP S 
mSl 2 



h(l, si) 



W FF S (L 2 
mSl 2 \ n2 



(Bl) 



L 

u c 



We can replace spectral frequency SI with uj gw since 
h(St) ~ hoS(Sl ± w gw ). Then 



X(Sl) ~Lh{Sl) 



mSl 2 



£fp 

mc 2 



1 



A 



I) Lh(Q), 



where £fp = WppSL is the full energy of light wave in- 
side the cavity. We can neglect the third term in formula 
dEJ completely since £ FP /mc 2 w 10~ 17 (£ FP = 20 J 
and m = 40 kg for Advanced LIGO project), Ao/A gw ~ 
10~ 13 ~ 10- 10 and A /L « 10" 9 . 

Now looking back, we see that the term we have just 
omitted corresponds to Fq^ — the statement required 
to be proved. 



APPENDIX C: CALCULATION OF THE FORCE 
OF LIGHT PRESSURE IN AN EXPLICIT FORM 

In this Appendix we calculate explicitly the force of 
light pressure 



Q 

FJX,X,t) = — 



1 fdA 
dt 



OA 

dx 



(i) 



x=L+X(t) 



which is the right side of Eq. ([23]) (for simplicity we use 
notation F x instead of F^). 



1. General formula 

First we calculate the force of light pressure in time 
domain. Summing up the derivatives of A + (x,t) and 
A_(x,t) (see formula (jllap ). raising the sum to the sec- 
ond power, neglecting the terms containing e 2luJat and 



keeping only the first order in h terms, we obtain: 



OA 

dx 



L+X{t) 



= 2ik Aj, 

+A +0 
-A-o 
-A-o 

+C.C. 



Substituting formula A+o 
tain: 



( da\ 



\ dx 
da* 



dx 
da*_ 

dx 
da\ 

dx 



+ ik Q a*_ ) e 2ikoL 



+ ikoa*_ 

- ik Q a{] e - 2lkaL 



dA 

dx 



L+X(t) 



^ik A +0 
-Mk A- 



-A-qc 2lU0T here we ob- 
dal » 

ikoa'L 



dx 

da*_ 
dx 



c.c. 
(CI) 



Similar calculations lead to 



1 fdA 

at 



L+X{t) 



This is not surprising since E Z (L + X) = —(X /c)H y (L + 
X) (see formula (fl4|) ) and being squared leads to X 2 . 
Therefore, E 2 equals to zero in our linear approximation. 

Substituting formula (jllbp into (|C1[) we obtain the fi- 
nal formula for _F r : 



F x (Sl) 



F x {X,X,t) = I F x {Sl)e 
Sk 



2tt 



A* +0 a + (u] + SI) 
+A +Q a* + (ujo - Si) 
+A*_ a_{Lj a + Si) 

+A_ a*L{uJ - SI) 



w + 


si 


c 






si 


c 




UJO + 


si 


c 




UJQ - 


si 



-iat dn 

2tt' 



JUt 



(C2) 



Now we substitute explicit formulas (|17a|) and (|17b[) here 
and split F x (SI) into several parts according to their phys- 
ical sense: 

f x (si) = F Ruct (si) + F pm (si) + (SI). 

Note that we have not taken into account the terms 
proportional to A±og±(x, t) while calculating the force of 
light pressure. It is straightforward to verify, using sim- 
ilar calculations, that the corresponding correction can 
be estimated as SF x (Sl) w SW FP (SlT) 2 h(Sl), resulting in 
inaccuracy SX(Sl) « (£ F p /mc 2 )Lh(Sl) « 10~ 17 Lh(Sl). 



16 



2. Fluctuating force 

Here we obtain formula for the fluctuating part 
-Ffluct(^) of the radiation pressure arising due to fluc- 
tuations a; n (x, t) of the input light wave: 



Ffluct(ft) 

Sk 



2- 



ain(wo + fi) 

+a* n (w - fi) 
-a in (w + fi) 
-a* n (w - fi) 



1 — i? e 2i(w +O)T c 

X — Re~ 2i(wa— H)t c 
j4 ^ oTe 2i( Wo +n)r e - l Or Wf) + ^ 
X _ _^ e 2i(a) +Sl)r c 



1 _ J£g-2i(aJ -O)T 

Using formulas (|16ap and (|16b|) we obtain: 



STeo 

7TC 



T 2 4* 



T 2 An (w - n)e'" T _ 

1 _ Re^or l _ i?e -2 i (o;o-n)r a in\ UJ "J 



3. Ponderomotive force 

Here we obtain the formula for the coefficient /C(fi) 
and corresponding ponderomotive force F pm (Sl) = 
-/C(Q)X(fi) (see formulas (g4j) and Ponderomo- 
tive force is the sum of all the terms in formula (IC2|) that 
include phase shift (|18|) arising due to mirror's displace- 
ment X (see formulas (|17a[) and (jTYbJ) ) : 



pm 



*1 
SW FP 

" 2fc 



Re 2iu ° T e inT 



w 



i<5* mir (fi) 



FP 



2fco 



Re -2iu r e iQr 
1 _ Re -2i(ui -Q,)T ' i _ Re -2i(ui -Q)7 



W 



Substituting formula (|18p here we obtain the following 
formula for /C(fi): 



-i5W 



FP 



1 + /£ e 2i ( ClJo + f2 ) T 



w 



1 — R e 2i ( UJ a+ n ) T c 



I _ Re -2i( Uo -n)r c a - i] 

It will be convenient for the future prospects to divide 
fC(Q) into two parts: £(fi) = K(O) - 2iOT(0) using 



the (wo ± f2)/c multipliers. The first part K(fl) is pro- 
portional to wo I c and is called the coefficient of optical 
rigidity: 



K(0) 



Ak SW FP Re 2lUT sin2w T 
1 - 2Re 2inT cos 2w r + R 2 e ilCtT ' 



The second part — 2if2r(f2) is proportional to f2/c where 
r(f2) is called the coefficient of radiative friction: 



r(0) = 



5Wi 



FP 



l-R z e 



c 1 - 2Re 2mT cos 2w r + R 2 e uo - T ' 



Therefore, we secure full analogy with linear harmonic 
oscillator's mechanical impedance. Obviously, |2f2r| ~ 

Wvo)\K\- 



GW+EMW correction 



Here we obtain formula ([26]) for the 'GW+EMW cor- 
rection to the force of light pressure. It arises due to di- 
rect interaction of gravitational wave with the light wave 
inside the cavity and becomes significant when length of 
the cavity becomes comparable with gravitational wave- 
length. Therefore, this part of (|C2[) is composed of the 
terms that include phase shift (fT9)) arising due to direct 
coupling 'GW+EMW: 



F, 



gw+emw(^) 



FP 



2fc 



Re 2iuJ T e mT 



1 — J{ e 2i( UJ o+Q) T 1 — ^ e 2i(LJo+n)T 



2 ^gw+cmw(^) 



SWi 



FP 



i?e- 2i " 0T e inT 

I _ ^ e -2i(w -n)T J _ Re -2i(uj -Q)r 

u>o — Q . 
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gw+cmw 



(-0). 



Substituting formula (|20|) here we obtain: 



gw+omw(^) — — — Lh(Q.) ( 1 — 



sin fh 



1 



fir 



FP 



1 — /j; e 2i ( Wo+n ) T c 
1 + j R e - 2i ^°-^ T w 



1 _ Re -2i{ui -Q,)T 



Taking into account formula (|C3p we finally obtain the 
following formula for the 'GW+EMW correction: 

F gw+ew (n) = ix:(n) (1 - SJ ^) Lh(ci), 

which coincides with (|26p . 
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